Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies.

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

The chickpea WIP2 gene underlying a major QTL contributes to lateral root development.

Lateral roots are a major component of root system architecture, and lateral root count (LRC) positively contributes to yield under drought in chickpea. To understand the genetic regulation of LRC, a biparental mapping population derived from two chickpea accessions having contrasting LRCs was genotyped by sequencing, and phenotyped to map four major quantitative trait loci (QTLs) contributing to 13-32% of the LRC trait variation. A single- nucleotide polymorphism tightly linked to the locus contributing to highest trait variation was located on the coding region of a gene (CaWIP2), orthologous to NO TRANSMITTING TRACT/WIP domain protein 2 (NTT/WIP2) gene of Arabidopsis thaliana. A polymorphic simple sequence repeat (SSR) in the CaWIP2 promoter showed differentiation between low versus high LRC parents and mapping individuals, suggesting its utility for marker-assisted selection. CaWIP2 promoter showed strong expression in chickpea apical root meristem and lateral root primordia. Expression of CaWIP2 under its native promoter in the Arabidopsis wip2wip4wip5 mutant rescued its rootless phenotype to produce more lateral roots than the wild-type plants, and led to formation of amyloplasts in the columella. CaWIP2 expression also induced the expression of genes that regulate lateral root emergence. Our study identified a gene-based marker for LRC which will be useful for developing drought-tolerant, high-yielding chickpea varieties.

Assessment of diversity in anti-nutrient profile, resistant starch, minerals and carbohydrate components in different ricebean (Vigna umbellata) accessions.

Ricebean accessions (n = 38) cultivated in India were evaluated for their comprehensive nutrient, anti-nutrients and mineral composition. Protein and total dietary fibre ranged between 23.23 and 27.33 and 12.27 to 16.69 g/100 g, respectively. Among the oligosaccharides, verbascose was not detected, however, raffinose and stachyose ranged between 47 and 186 and 117 to 5765 mg/100 g, respectively. Among the free sugars, sucrose was found dominating (up to 370 mg/100 g). Resistant starch (4.13 to 8.62 %), iron (3.49 to 7.46 mg/100 g), zinc (1.90 to 3.72 mg/100 g) and selenium (0.28 to 4.48 µg/100 g) varied significantly (p < 0.05) among ricebean samples. Phytic acid, saponin, trypsin inhibitor and oxalate analysed in ricebean accessions ranged between 303 and 760 mg/100 g, 19 to 46 mg/g, 309 to 1076 mg/100 g and 219 to 431 mg/100 g, respectively. Multivariate analysis using hierarchical clustering analysis (HCA), and principal component analysis (PCA) was employed to decipher the diversity of nutrients and anti-nutrients across the ricebean accessions. Based on HCA, dendrogram-1 (nutrients) and dendrogram-2 (minerals, anti-nutrients) were produced, having four clusters in each. In the dendrogram-1 and 2, the largest cluster had (n = 21) and (n = 15) accessions, respectively. The PCA analyse the uncorrelated set of variables (principal components) and it condenses a large set of data variables. Based on the eigenvalue >1, a total of eight PCs were formed contributing total variance of 78.8 %. The factor loading contribution in the PC1 and PC2 were from iron, fructose, glucose, raffinose and total dietary fibre, selenium (Se) and protein, respectively.

Response of Physiological, Reproductive Function and Yield Traits in Cultivated Chickpea (Cicer arietinum L.) Under Heat Stress.

Under global climate change, high-temperature stress is becoming a major threat to crop yields, adversely affecting plant growth, and ultimately resulting in significant yield losses in various crops, including chickpea. Thus, identifying crop genotypes with increased heat stress (HS) tolerance is becoming a priority for chickpea research. Here, we assessed the response of seven physiological traits and four yield and yield-related traits in 39 chickpea genotypes grown in normal-sown and late-sown environments [to expose plants to HS (>32/20°C) at the reproductive stage] for two consecutive years (2017-2018 and 2018-2019). Significant genetic variability for the tested traits occurred under normal and HS conditions in both years. Based on the tested physiological parameters and yield-related traits, GNG2171, GNG1969, GNG1488, PantG186, CSJ515, RSG888, RSG945, RVG202, and GNG469 were identified as promising genotypes under HS. Further, ten heat-tolerant and ten heat-sensitive lines from the set of 39 genotypes were validated for their heat tolerance (32/20°C from flowering to maturity) in a controlled environment of a growth chamber. Of the ten heat-tolerant genotypes, GNG1969, GNG1488, PantG186, RSG888, CSJ315, and GNG1499 exhibited high heat tolerance evidenced by small reductions in pollen viability, pollen germination, and pod set %, high seed yield plant-1 and less damage to membranes, photosynthetic ability, leaf water status, and oxidative processes. In growth chamber, chlorophyll, photosynthetic efficiency, pollen germination, and pollen viability correlated strongly with yield traits. Thus, GNG1969, GNG1488, PantG186, RSG888, CSJ315, and GNG1499 genotypes could be used as candidate donors for transferring heat tolerance traits to high-yielding heat-sensitive varieties to develop heat-resilient chickpea cultivars.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata).

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development-related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development-related traits in ricebeans. A total of 51 genes encoding SCF TIR1/AFB , Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean-derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

Breeding, Genetics, and Genomics Approaches for Improving Fusarium Wilt Resistance in Major Grain Legumes.

Fusarium wilt (FW) disease is the key constraint to grain legume production worldwide. The projected climate change is likely to exacerbate the current scenario. Of the various plant protection measures, genetic improvement of the disease resistance of crop cultivars remains the most economic, straightforward and environmental-friendly option to mitigate the risk. We begin with a brief recap of the classical genetic efforts that provided first insights into the genetic determinants controlling plant response to different races of FW pathogen in grain legumes. Subsequent technological breakthroughs like sequencing technologies have enhanced our understanding of the genetic basis of both plant resistance and pathogenicity. We present noteworthy examples of targeted improvement of plant resistance using genomics-assisted approaches. In parallel, modern functional genomic tools like RNA-seq are playing a greater role in illuminating the various aspects of plant-pathogen interaction. Further, proteomics and metabolomics have also been leveraged in recent years to reveal molecular players and various signaling pathways and complex networks participating in host-pathogen interaction. Finally, we present a perspective on the challenges and limitations of high-throughput phenotyping and emerging breeding approaches to expeditiously develop FW-resistant cultivars under the changing climate.

Salinity stress response and 'omics' approaches for improving salinity stress tolerance in major grain legumes.

KEY MESSAGE: Sustaining yield gains of grain legume crops under growing salt-stressed conditions demands a thorough understanding of plant salinity response and more efficient breeding techniques that effectively integrate modern omics knowledge. Grain legume crops are important to global food security being an affordable source of dietary protein and essential mineral nutrients to human population, especially in the developing countries. The global productivity of grain legume crops is severely challenged by the salinity stress particularly in the face of changing climates coupled with injudicious use of irrigation water and improper agricultural land management. Plants adapt to sustain under salinity-challenged conditions through evoking complex molecular mechanisms. Elucidating the underlying complex mechanisms remains pivotal to our knowledge about plant salinity response. Improving salinity tolerance of plants demand enriching cultivated gene pool of grain legume crops through capitalizing on 'adaptive traits' that contribute to salinity stress tolerance. Here, we review the current progress in understanding the genetic makeup of salinity tolerance and highlight the role of germplasm resources and omics advances in improving salt tolerance of grain legumes. In parallel, scope of next generation phenotyping platforms that efficiently bridge the phenotyping-genotyping gap and latest research advances including epigenetics is also discussed in context to salt stress tolerance. Breeding salt-tolerant cultivars of grain legumes will require an integrated "omics-assisted" approach enabling accelerated improvement of salt-tolerance traits in crop breeding programs.

mQTL-seq and classical mapping implicates the role of an AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED (AHL) family gene in Ascochyta blight resistance of chickpea.

Ascochyta blight (AB) caused by the fungal pathogen Ascochyta rabiei is a serious foliar disease of chickpea (Cicer arietinum L.). Despite many genetic studies on chickpea-Ascochyta interaction, genome-wide scan of chickpea for the identification of AB-associated quantitative trait loci (QTLs) and their gene(s) has not been accomplished. To elucidate narrow QTLs for AB resistance, here, we report the use of multiple QTL-sequencing approach on 2 sets of extreme AB phenotype bulks derived from Cicer intraspecific and interspecific crosses. Two major QTLs, qABR4.1 and qABR4.2, and a minor QTL, qABR4.3, were identified on assembled chickpea pseudomolecule 4. We narrowed qABR4.1 to a "robust region" at 4.568-4.618 Mb through mapping on a larger intraspecific cross-derived population and comparative analysis. Among 4 genes, the CaAHL18 gene showed higher expression under Ascochyta stress in AB resistant parent suggesting that it is the candidate gene under "robust qABR4.1." Dual-luciferase assay with CaAHL18 polymorphic cis-regulatory sequences showed that allelic variation is associated with higher expression. Thus, our findings on chickpea-Ascochyta interaction have narrowed down AB resistance associated QTLs on chickpea physical map. The narrowed QTLs and gene-associated markers will help in biotechnological and breeding programs for chickpea improvement.

Population structure and association analysis of heat stress relevant traits in chickpea (Cicer arietinum L.).

Understanding genetic diversity and population structure is prerequisite to broaden the cultivated base of any crop. In the current investigation, we report discovery of a total of 319 alleles by assaying 81 SSRs on 71 chickpea genotypes. The cluster analysis based on Jaccard coefficient and unweighted neighbor joining algorithm categorized all genotypes into two major clusters. Cultivars grown within the same agro-climatic zones were clustered together, whereas the remaining genotypes particularly advanced breeding lines and accessions assigned to another cluster. Population structure analysis separated the entire collection into two subpopulations (K = 2) and the clustering pattern remained in close agreement with those of distance-based methods. Importantly, we also discovered marker trait association for membrane stability index (MSI) and leaf chlorophyll content measured as SPAD chlorophyll meter reading (SCMR), the two important physiological parameters indicative of heat stress (HS) tolerance in chickpea. Association analysis using both general linear and mixed linear models of the mean phenotypic data of traits recorded in 2016 and 2017 uncovered significant association of NCPGR206 and H2L102 with the MSI trait. Likewise, SSR markers GA9, TR31 and TA113 exhibited significant association with SCMR trait. The genomic regions putatively linked with two traits may be investigated in greater detail to further improve knowledge about the genetic architecture of HS tolerance in chickpea.

A repeat length variation in myo-inositol monophosphatase gene contributes to seed size trait in chickpea.

Chickpea (Cicer arietinum L.) is the third most important food legume crop. Seed size is the most economically important trait for chickpea. To understand the genetic regulation of seed size in chickpea, the present study established a three-way association of CT repeat length variation of a simple sequence repeat (SSR) in myo-inositol monophosphatase gene (CaIMP) with seed weight and phytic acid content by large scale validation and genotyping in a set of genetically diverse germplasm accessions and two reciprocal intra-specific mapping populations. Germplasms and mapping individuals with CT repeat-length expansion in the 5' untranslated region of CaIMP exhibited a pronounced increase in CaIMP protein level, enzymatic activity, seed-phytate content and seed weight. A chickpea transient expression system demonstrated this repeat-length variation influenced the translation of CaIMP mRNA, apparently by facilitating translation initiation. Our analyses proposed that the SSR marker derived from 5' UTR of a CaIMP gene is a promising candidate for selection of seed size/weight for agronomic trait improvement of chickpea.

Genome-wide investigation and expression analyses of WD40 protein family in the model plant foxtail millet (Setaria italica L.).

WD40 proteins play a crucial role in diverse protein-protein interactions by acting as scaffolding molecules and thus assisting in the proper activity of proteins. Hence, systematic characterization and expression profiling of these WD40 genes in foxtail millet would enable us to understand the networks of WD40 proteins and their biological processes and gene functions. In the present study, a genome-wide survey was conducted and 225 potential WD40 genes were identified. Phylogenetic analysis categorized the WD40 proteins into 5 distinct sub-families (I-V). Gene Ontology annotation revealed the biological roles of the WD40 proteins along with its cellular components and molecular functions. In silico comparative mapping with sorghum, maize and rice demonstrated the orthologous relationships and chromosomal rearrangements including duplication, inversion and deletion of WD40 genes. Estimation of synonymous and non-synonymous substitution rates revealed its evolutionary significance in terms of gene-duplication and divergence. Expression profiling against abiotic stresses provided novel insights into specific and/or overlapping expression patterns of SiWD40 genes. Homology modeling enabled three-dimensional structure prediction was performed to understand the molecular functions of WD40 proteins. Although, recent findings had shown the importance of WD40 domains in acting as hubs for cellular networks during many biological processes, it has invited a lesser research attention unlike other common domains. Being a most promiscuous interactors, WD40 domains are versatile in mediating critical cellular functions and hence this genome-wide study especially in the model crop foxtail millet would serve as a blue-print for functional characterization of WD40s in millets and bioenergy grass species. In addition, the present analyses would also assist the research community in choosing the candidate WD40s for comprehensive studies towards crop improvement of millets and biofuel grasses.

Population structure and association mapping of yield contributing agronomic traits in foxtail millet.

Association analyses accounting for population structure and relative kinship identified eight SSR markers ( p < 0.01) showing significant association ( R (2) = 18 %) with nine agronomic traits in foxtail millet. Association mapping is an efficient tool for identifying genes regulating complex traits. Although association mapping using genomic simple sequence repeat (SSR) markers has been successfully demonstrated in many agronomically important crops, very few reports are available on marker-trait association analysis in foxtail millet. In the present study, 184 foxtail millet accessions from diverse geographical locations were genotyped using 50 SSR markers representing the nine chromosomes of foxtail millet. The genetic diversity within these accessions was examined using a genetic distance-based and a general model-based clustering method. The model-based analysis using 50 SSR markers identified an underlying population structure comprising five sub-populations which corresponded well with distance-based groupings. The phenotyping of plants was carried out in the field for three consecutive years for 20 yield contributing agronomic traits. The linkage disequilibrium analysis considering population structure and relative kinship identified eight SSR markers (p < 0.01) on different chromosomes showing significant association (R (2) = 18 %) with nine agronomic traits. Four of these markers were associated with multiple traits. The integration of genetic and physical map information of eight SSR markers with their functional annotation revealed strong association of two markers encoding for phospholipid acyltransferase and ubiquitin carboxyl-terminal hydrolase located on the same chromosome (5) with flag leaf width and grain yield, respectively. Our findings on association mapping is the first report on Indian foxtail millet germplasm and this could be effectively applied in foxtail millet breeding to further uncover marker-trait associations with a large number of markers.

Development of 5123 intron-length polymorphic markers for large-scale genotyping applications in foxtail millet.

Generating genomic resources in terms of molecular markers is imperative in molecular breeding for crop improvement. Though development and application of microsatellite markers in large-scale was reported in the model crop foxtail millet, no such large-scale study was conducted for intron-length polymorphic (ILP) markers. Considering this, we developed 5123 ILP markers, of which 4049 were physically mapped onto 9 chromosomes of foxtail millet. BLAST analysis of 5123 expressed sequence tags (ESTs) suggested the function for ∼71.5% ESTs and grouped them into 5 different functional categories. About 440 selected primer pairs representing the foxtail millet genome and the different functional groups showed high-level of cross-genera amplification at an average of ∼85% in eight millets and five non-millet species. The efficacy of the ILP markers for distinguishing the foxtail millet is demonstrated by observed heterozygosity (0.20) and Nei's average gene diversity (0.22). In silico comparative mapping of physically mapped ILP markers demonstrated substantial percentage of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (∼50%), maize (∼46%), rice (∼21%) and Brachypodium (∼21%) chromosomes. Hence, for the first time, we developed large-scale ILP markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.

Genetic diversity and population structure of rice landraces from Eastern and North Eastern States of India.

BACKGROUND: Adaptations to different habitats across the globe and consequent genetic variation within rice have resulted in more than 120,000 diverse accessions including landraces, which are vital genetic resources for agronomic and quality traits. In India the rice landraces of the states West Bengal, Assam, Mizoram, Manipur and Nagaland are worthy candidates for genetic assessment. Keeping the above in view, the present study was conducted with the aim to (i) calculate the genetic distances among the accessions of 83 landraces collected from these states along with 8 check accessions (total 91 accessions) using 23 previously mapped SSR markers and (ii) examine the population structure among the accessions using model-based clustering approach. RESULTS: Among the 91 accessions, 182 alleles were identified which included 51 rare and 27 null alleles. The average PIC value was 0.7467/marker. The non-aromatic landraces from West Bengal was most diverse with 154 alleles and an average PIC value of 0.8005/marker, followed by the aromatic landraces from West Bengal with 118 alleles and an average PIC value of 0.6524/marker, while the landraces from North East ranked third with 113 alleles and an average PIC value of 0.5745/marker. In the dendrogram distinct clusters consisting of predominantly aromatic landraces and predominantly North East Indian landraces were observed. The non-aromatic landraces from West Bengal were interspersed within these two clusters. The accessions were moderately structured, showing four sub-populations (A-D) with an Fst value of 0.398, 0.364, 0.206 and 0.281, respectively. The assigned clustering of accessions was well in agreement in both distance-based and model-based approaches. CONCLUSIONS: Each of the accessions could be identified unequivocally by the SSR profiles. Genetically the non aromatic landraces from West Bengal were most diverse followed by the aromatic landraces from the same state. The North Eastern accessions ranked third. Further, grouping of accessions based on their agronomic traits may serve as a resource for future studies, leading to the improvement of rice. Moreover in-situ preservation of the landraces is also a means of protection of biodiversity and cultural heritage.

Development of eSSR-Markers in Setaria italica and Their Applicability in Studying Genetic Diversity, Cross-Transferability and Comparative Mapping in Millet and Non-Millet Species.

Foxtail millet (Setariaitalica L.) is a tractable experimental model crop for studying functional genomics of millets and bioenergy grasses. But the limited availability of genomic resources, particularly expressed sequence-based genic markers is significantly impeding its genetic improvement. Considering this, we attempted to develop EST-derived-SSR (eSSR) markers and utilize them in germplasm characterization, cross-genera transferability and in silico comparative mapping. From 66,027 foxtail millet EST sequences 24,828 non-redundant ESTs were deduced, representing ~16 Mb, which revealed 534 (~2%) eSSRs in 495 SSR containing ESTs at a frequency of 1/30 kb. A total of 447 pp were successfully designed, of which 327 were mapped physically onto nine chromosomes. About 106 selected primer pairs representing the foxtail millet genome showed high-level of cross-genera amplification at an average of ~88% in eight millets and four non-millet species. Broad range of genetic diversity (0.02-0.65) obtained in constructed phylogenetic tree using 40 eSSR markers demonstrated its utility in germplasm characterizations and phylogenetics. Comparative mapping of physically mapped eSSR markers showed considerable proportion of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (~68%), maize (~61%) and rice (~42%) chromosomes. Synteny analysis of eSSRs of foxtail millet, rice, maize and sorghum suggested the nested chromosome fusion frequently observed in grass genomes. Thus, for the first time we had generated large-scale eSSR markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.

Comprehensive genome-wide survey, genomic constitution and expression profiling of the NAC transcription factor family in foxtail millet (Setaria italica L.).

The NAC proteins represent a major plant-specific transcription factor family that has established enormously diverse roles in various plant processes. Aided by the availability of complete genomes, several members of this family have been identified in Arabidopsis, rice, soybean and poplar. However, no comprehensive investigation has been presented for the recently sequenced, naturally stress tolerant crop, Setaria italica (foxtail millet) that is famed as a model crop for bioenergy research. In this study, we identified 147 putative NAC domain-encoding genes from foxtail millet by systematic sequence analysis and physically mapped them onto nine chromosomes. Genomic organization suggested that inter-chromosomal duplications may have been responsible for expansion of this gene family in foxtail millet. Phylogenetically, they were arranged into 11 distinct sub-families (I-XI), with duplicated genes fitting into one cluster and possessing conserved motif compositions. Comparative mapping with other grass species revealed some orthologous relationships and chromosomal rearrangements including duplication, inversion and deletion of genes. The evolutionary significance as duplication and divergence of NAC genes based on their amino acid substitution rates was understood. Expression profiling against various stresses and phytohormones provides novel insights into specific and/or overlapping expression patterns of SiNAC genes, which may be responsible for functional divergence among individual members in this crop. Further, we performed structure modeling and molecular simulation of a stress-responsive protein, SiNAC128, proffering an initial framework for understanding its molecular function. Taken together, this genome-wide identification and expression profiling unlocks new avenues for systematic functional analysis of novel NAC gene family candidates which may be applied for improvising stress adaption in plants.

Genome-wide development and use of microsatellite markers for large-scale genotyping applications in foxtail millet [Setaria italica (L.)].

The availability of well-validated informative co-dominant microsatellite markers and saturated genetic linkage map has been limited in foxtail millet (Setaria italica L.). In view of this, we conducted a genome-wide analysis and identified 28 342 microsatellite repeat-motifs spanning 405.3 Mb of foxtail millet genome. The trinucleotide repeats (∼48%) was prevalent when compared with dinucleotide repeats (∼46%). Of the 28 342 microsatellites, 21 294 (∼75%) primer pairs were successfully designed, and a total of 15 573 markers were physically mapped on 9 chromosomes of foxtail millet. About 159 markers were validated successfully in 8 accessions of Setaria sp. with ∼67% polymorphic potential. The high percentage (89.3%) of cross-genera transferability across millet and non-millet species with higher transferability percentage in bioenergy grasses (∼79%, Switchgrass and ∼93%, Pearl millet) signifies their importance in studying the bioenergy grasses. In silico comparative mapping of 15 573 foxtail millet microsatellite markers against the mapping data of sorghum (16.9%), maize (14.5%) and rice (6.4%) indicated syntenic relationships among the chromosomes of foxtail millet and target species. The results, thus, demonstrate the immense applicability of developed microsatellite markers in germplasm characterization, phylogenetics, construction of genetic linkage map for gene/quantitative trait loci discovery, comparative mapping in foxtail millet, including other millets and bioenergy grass species.

Genic non-coding microsatellites in the rice genome: characterization, marker design and use in assessing genetic and evolutionary relationships among domesticated groups.

BACKGROUND: Completely sequenced plant genomes provide scope for designing a large number of microsatellite markers, which are useful in various aspects of crop breeding and genetic analysis. With the objective of developing genic but non-coding microsatellite (GNMS) markers for the rice (Oryza sativa L.) genome, we characterized the frequency and relative distribution of microsatellite repeat-motifs in 18,935 predicted protein coding genes including 14,308 putative promoter sequences. RESULTS: We identified 19,555 perfect GNMS repeats with densities ranging from 306.7/Mb in chromosome 1 to 450/Mb in chromosome 12 with an average of 357.5 GNMS per Mb. The average microsatellite density was maximum in the 5' untranslated regions (UTRs) followed by those in introns, promoters, 3'UTRs and minimum in the coding sequences (CDS). Primers were designed for 17,966 (92%) GNMS repeats, including 4,288 (94%) hypervariable class I types, which were bin-mapped on the rice genome. The GNMS markers were most polymorphic in the intronic region (73.3%) followed by markers in the promoter region (53.3%) and least in the CDS (26.6%). The robust polymerase chain reaction (PCR) amplification efficiency and high polymorphic potential of GNMS markers over genic coding and random genomic microsatellite markers suggest their immediate use in efficient genotyping applications in rice. A set of these markers could assess genetic diversity and establish phylogenetic relationships among domesticated rice cultivar groups. We also demonstrated the usefulness of orthologous and paralogous conserved non-coding microsatellite (CNMS) markers, identified in the putative rice promoter sequences, for comparative physical mapping and understanding of evolutionary and gene regulatory complexities among rice and other members of the grass family. The divergence between long-grained aromatics and subspecies japonica was estimated to be more recent (0.004 Mya) compared to short-grained aromatics from japonica (0.006 Mya) and long-grained aromatics from subspecies indica (0.014 Mya). CONCLUSION: Our analyses showed that GNMS markers with their high polymorphic potential would be preferred candidate functional markers in various marker-based applications in rice genetics, genomics and breeding. The CNMS markers provided encouraging implications for their use in comparative genome mapping and understanding of evolutionary complexities in rice and other members of grass family.